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ON  THE  ACCELERATION  OF  MAGNETOGASDYNAMIC  CHANNEL 
FLOWS  THROUGH  THE  SONIC  VELOCITY 

by  F.  D.  Hains 


The  possibility  of  passage  through  the  sonic 
velocity  is  examined  for  the  magnetogasdynamic  flow  through 
a  diverging,  straight-walled  channel.  In  addition  to  the 
usual  assumptions  of  continuum  MGD,  the  electrical  conduc¬ 
tivity  is  assumed  to  be  small,  but  the  effect  of  Hall 
current  is  included.  A  number  of  different  cases  arise, 
depending  on  whether  or  not  one  or  more  parameters  vanish 
or  are  small. 

Numerical  results  for  the  case  where  the  Hall 
current  vanishes,  and  the  problem  becomes  analogous  to  the 
gas dynamic  source,  indicate  the  sonic  point  can  be  placed 
anywhere  along  the  channel  by  proper  selection  of  the 
electric  and  magnetic  field  strengths  and  boundary  condi¬ 


tions 


I.  INTRODUCTION 


This  report  deals  with  the  flow  of  compressible,  electri¬ 
cally  conducting  gases  of  sufficient  density  to  satisfy  the  continuum 
assumptions  of  gasdynamics.  An  attempt  will  be  made  to  obtain  exact 
solutions,  either  closed  form  or  numerical,  of  the  magnetogasdynamic 
flow  equations*  Since  the  MGD  equations  reduce  to  those  of  ordinary 
gas  dynamics  when  the  electrical  conductivity  vanishes,  it  is  logi¬ 
cal  to  seek  solutions  which  will  reduce  to  known  exact  solutions  of 
gasdynamics.  Two  exact  solutions  of  gaBdynamics  ax;e  the  source- 
sink  and  source-vortex.  This  report  will  consider  the  MGD' source 
and  will  emphasize  those  solutions  which  pass  through  the  sonic 
velocity. 

In  Reference  1,  three  simple  solutions  for  the  MGD  source 

were. presented.  The  configuration  studied  is  shown  in  Figure  1.  A 
,  *  \ 
source  flow  with  radial  velocity  u  interacts  with  crossed  electric 

and  magnetic  field.  The  electric  field  is  uniform  in  the  z  direc¬ 
tion,  while  the  magnetic  field  is  produced  by  a  wire  along  the  z 
axis.  The  three  cases  considered  are  summarized  in  Table  I. 


TABLE  I 


Case 

R 

m 

.  K 

E 

P 

Comments 

A 

,.P 

0 

0 

0 

Finite 

Closed  form  solution 

B 

0 

0 

Finite 

0 

Solution  sketched 

C 

0 

0 

Finite 

Finite 

One  particular  solution 

Closed  form  solutions  were  obtained  in  Case  A,  where  the  electric 


i 
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field  vanishes  and  the  flow  is  isoenergetic.  The  velocity  distri¬ 
butions  for  various  values  of  p  are  shown  in  Figure  3»  Only  the 
solution  for  p  =  0.6  passes  through  the  sonic  velocity  u  *  0.5* 

All  other  solutions  terminate  with  sonic  velocity  at  a  limit  circle. 
Some  of  these  solutions  lie  inside  a  limit  circle  while  others  lie 
outside  a  limit  circle.  The  solutions  for  Case  B  are  sketched  in 
Figure  4.  It  appears  that  all  solutions  for  this  case  are  bounded 
in  finite  regions  formed  by  two  limit  circles.  In  Case  C,  a  parti¬ 
cular  solution  was  obtained  which  proceeds  at  constant  Mach  number. 

The  cases  discussed  in  this  report  are  shown  in  Table  II. 
Because  of  the  non-linearity  of  the  equations,  numerical  integra¬ 
tion  is  usually  required  for  these  cases. 


TABLE  II 


Case 

R 

m 

K 

E  . 

P 

Comments 

1 

0 

Finite 

Finite 

Finite 

2 

0 

Finite 

Finite 

Finite 

One  particular  solution 

3 

0 

'  0 

Finite 

Finite 

w  ^  0 

4 

0 

0 

Finite 

Finite 

ii 

o 

5 

Finite 

0 

0 

Finite 

6 

Finite 

0 

0 

0 

Solve  for  Bq 

II.  BASIC  EQUATIONS 
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The  basic  MGD  equations  in  vector  form  are 


continuity 

V.pq  s  0 

(1) 

momentum 

p(q*V)q  +  Vp  =  j*  x  B 

(2) 

energy 

p(q*V)H  =  E  •  j* 

(3) 

Ohm ' s  law 

j  =  o[e  +  q  x  B  -  jkjf] 

e 

(4) 

Maxwell ' s 
equations 

V  .  j  =  0 

(5) 

j  =  V  X  B 

(6) 

V  .  B  =  0 

(7) 

<n 

X 

Mi 

II 

o 

(8) 

V  »  E  =  0 

(9) 

gas  law 

p  =  pRT 

(10) 

Equations  (l)  -  (9)  can  be  expressed  in  the  following  scalar  form  in 
the  r-9-z  coordinate  system  of  Figure  1. 


r  a7(rpu)  +r  dVpv)  +  a7(pw)  =  0 


Tu  £  +  w  »«1+  '*£  = 

|_  dr  r  r  dzj  dr 

[dv  v  dv  uv  dv  dp 

u  d7+7d3+T  +  w  dlJ+Ss 

[" 


JeV  ^Be 


^zBr  ’  JrB. 


dw  v  dw  dw 
dr  +  r  T5  +  w  dz 


+  4*  = 


dz  "  jrB9  “  j0Br 


(11) 

(12) 

(15) 


(14) 
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dH  ov  dH  dH 

pu  d7+^"  d£  +  pw  d7  =  Erjr  +  E0J0  +  Vz 

(15) 

ir  =  -  JzBe)] 

(16) 

le  *  0 [_E0  +  wBr  ~ uBz  ” n^T^zBr  ’  W] 

(17) 

iz  =  «[BI*ttB8-TBr-^j<jrB8  -  JeBr)] 

(18) 

id,..  i  aj‘e  diz  A 
;a?(rjr)+71T  +  IT  =  0 

(19) 

1/9Br  drBgV 
r \ d9  dz  y 

(20) 

,  /dB  dB  \ 

^0  -  F(ar  “  -a#) 

(21) 

•  _  i/drBe  .  dBr\ 

**z  “  r\  dr  d 0 J 

(22) 

,  a  ,  dB0  dB 

—  -r— (rB  )  +—  -r-jr-  +  — — 
r  dr  r  r  8  dz 


=  0 


(23) 


dEp  drE0 

IT  ’  — 58  0 


(24) 


dE  dE 

r  _  z 

dz  dr 


(25) 


drEg  dEr 

- - dT  =  0 


(26) 


.  drE  .  dEfl  dE 
1. _ r  1_  9  _ z 

r  dr  r  d  6  dz 


(27) 
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To  obtain  solutions  to  this  set  of  equations  without  further  simpli¬ 
fying  assumptions  would  be  difficult. 

For  the  source  shown  in  Figure  1,  we  can  assume 


Br  =  Bz  =  Er  =  Ee  a  aT  =  0 

E  =  E  (a  constant) 

z 


(28) 


Eqs.  (23)  -  (27)  are  automatically  satisfied  by  these  assumptions. 
In  addition,  Eq.  (13)  is  satisfied  by  v  =  0,  and  Eq.  (21)  leads  to 
jg  =  0.  The  remaining  equations  simplify  to 


— ’^-(rpu)  +rr— -( pw)  =  0 
r  dr  dz 


(29) 


du  du 

u  +  w  T” 

dr  dz 


+  ^  =  *3  B 


dr 


z"0 


(30) 


dw  dw 

U  -T—  +  W 

dr  dz 


*  4*  =  JB, 


dz 


(31) 


pu 


dH 

dr 


dH 

+  PW  dl  = 


(32). 


3t  *  °[-»b8  *  F7(WJ  l55) 

~(3rv] 

7 *  -k-  -  0  (»> 

jr  "  ~  r  A(rBe>  <56> 
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J2  *  r 


(37) 


Combining  Eq.  (33)  and  (34) 


j  =  - * - - 

“(0) 


OB 


9 


■  wBe  +  ITT (E+  uBe) 

e 


(38) 


and 


1  + 


■  oBgW 

E  +  uB«  +  — — 
0  n  e 
e 


(39) 


These  two  equations  are  used  to  eliminate  j  and  Jz  from  the  other 
equations. 

All  equations  up  to  this  point  are  dimensional,  although 
the  tilda  which  indicates  a  dimensional  quantity  has  been  omitted 
for  simplicity.  It  is  convenient  at  this  point  to  non-dimensionalize, 
using  the  following  definitions: 


/ 

r  =  r/r 


.  P  =  p/p„ 


u  =  u/um  p  =  p/p6 


=  z/r 


j  =  j  /au  B 
Jr  r  mo 


R  =  au  r 
m  mo 


T  =  T/T 


j  =  j  /au  B 
Jz  °z  mo 


N  = 


~  2 
aB  r 
o  o 

P  U 
s  m 


w  =  w/u 


m 


E  =  E/u  B 
m  o 


K 


aB 


n  e 
e 


B0  -  VBo 


H  =  SLO-1& 

/r-i  p„ 


m  =  m/rpu 
o  6  in 


(40) 


Many  of  these  are  the  same  used  in  Case  A  of  Reference  1  where  the 
flow  is  isoenergetic  (E  =  0)  so  that 


P  u 

Y  s  m 

p  =  1 2" 

6 


(41) 
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Although  E  is  generally  not  zero  in  this  report,  it  is  convenient 
to  retain  the  same  definitions. 

The  non-dimensional  form  of  the  basic  equations  is 


i-^rpu)  *  £<p»)  .  0 


du  Ju  So  -BB8[E.uBe.K.B|] 

p“  a? - ~XT 


1+K  B 


9 


dw  dw  dp 

pu  *r—  +  pw  —  +  a  = 
dr  dz  dz 


NBg|  -w+K(E+uB0)] 


2  2 
i+k2b| 


M  dH  2NE(E+uB0+KwB2) 


PU  "T—  +  pw  rr~  = 

dr  dz 


2  2 
1+K  Bi 


m 


i  d  ^Be("w+K[E+uBe]A  d  E+uBe+KwB 

r  dr\  , ,v2„2  I*  dz  1+K2B2 

9 


i*k2b| 


*  0 


is,.,  K»[»VKBe(E*»Be)] 

““  T“V  TiJ/s  y  —  -T 


r  dz'"  9 


i.k5b| 


14,  ,  Bm[E,uB9,K.BB]  • 


r  dr  8 


2  2 

1+K  B| 


where 


p  =  pT 


and 


H  =  T  +  u2 


2 

w 


(42) 

(43) 

(44) 

(45) 

(46) 

(47) 

(48) 

(49) 


(50) 


III.  CASE  1 


In  most  applications  of  MGD  channel  flow  the  electrical 
conductivity  is  small.  For  this  reason  this  case  will  be  restricted 
to  solutions  with  =  0.  Physically  this  means  the  magnetic  field 
will  interact  with  the  flow,  but  not  vice  versa.  The  basic  equations 
simplify  with  Eqs.  (47)  and  (48)  satisfied  by  B0  *  l/r,  and  the 
remaining  equations  become 

7  e7(rpu)  +  a!(Pw)  =  0  (51) 


PU  |SL  +  pw  +  a  =  -NC.Er  tu  +  Kw/r] 
dr  dz  dr  2  „2 

r  +K. 


(52) 


dw  dw  dp 

pu  -T—  +  pw  -r—  +  a  4*4 

dr  r  dz  <dz 


N[-w  +  K(E  + u/r) ] 


2  „2 

r  +  K 


pu 


dH  +  dH  i=  2NEr  [Eir. .+  u  t  Kw/r] 


dr 


dz 


r2  +IT 


(53) 

(54) 


In  this  form  the  equations  are  still  very  difficult  to  solve,  so  we 
d 

assume  rjy  =  0,  and  reduce  to  the  following  set  of  ordinary  differen¬ 
tial  equations 

rpu  =  m  (55) 


pu 


du  dp  -N[Er  + u  + Kw/r ] 

j t  _ 


dr 


dr 


2  2 
r  +  n 


(56) 


dw  N[-w  +  K(E  +  u/r)  ] 

dr  =  2  v2 

r  +  K 


dH  2NEr [Er  +  u  +  Kw/ r ] 

dr  =  2  v2 

r  +  K 


(57) 


(58) 
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Using  Eqs.  (49),  (50)  and  (55)  to  eliminate  p,  p,  and  H,  we  obtain 


du  _  F 
dr  -  G 


(59) 


where 


F  a  u< 


jp  r2u2  -  Er?  +  y^T  Kwr^u  “  2E2r^  -  2r2w2l+  (t2+K2)t| 

(60) 


G  =  r(r2+K2)  u2  -  t] 


(61) 


dw  _  JBj^  wr  +  K(  Er  +  u)  ] 


dr 


2  ,,2  > 


dT  _  2gEr2[Er  +  u  +  Kw/r]  _  du  -  dw 

dr  =  2  v2  ^  dr  ”  W  dr 

r  +  K 


(62) 

(65) 


For  given  boundary  conditions  and  specific  values  of  the  parameters 
Y,  P,  E  and  K,  these  equations  can  be  integrated  numerically  to 
yield  the  solution. 

Before  doing  this,  it  is  interesting  to  study  the  limit 
circle  where  F  /  0,  G  =  0.  Eq.  (59)  indicates  the  gradient  in  u  is 
infinite  so  the  solution  cannot  penetrate  the  limit  circle*  From  the 
definition  of  local  Mach  number  given  by 

» -  V(&)  ^ 

the  flow  must  be  supersonic  at  the  limit  circle.  As  in  Case  A  shown 
in  Figure  5,  we  expect  solutions  which  lie  inside  the  limit  circles, 
and  others  which  lie  outside  the  limit  circles.  The  type  of  results 
that  might  be  expected  are  sketched  in  Figure  5(a)  "  5(d). 
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In  this  report  we  seek  solutions  which  do  not  terminate  in 

limit  lines.  A  characteristic  point  (r  ,  u  ,  T  ,  w  )  is  one  which 

c  c  c  c 

satisfies  the  conditions  F  =  G  =  0.  Since  4“  is  indeterminate,  more 

dr 

than  one  integral  curve  may  pass  through  the  characteristic  point* 
The  singularity  condition  G  =  0  leads  to 


T  ■ 
c 


2 


u 


c 


(65) 


and  the  regularity  condition  F  =  0  leads  to  the  expression 


p|ru2  +  (2y-1+AW)U  +  (y-i)(i+w2)}  -Cl+\2)U2  a  0 


(66) 


In  the  last  equation  one  variable  has  been  eliminated  by  introducing 
the  quantities 


U  = 


Er 


W 


w 


Er 


rc 


(67) 


The  slopes  of  the  integral  curves  at  the  characteristic  point  are 


found  by  L' Hospital's  Rule. 

I du\  _  Lim  du_  _  Lim  dF/dr 
\dr  Jc  ~  r-rrc.  dr  ~  r-*rc  dG/dr  ' 


(68) 


Two  curves  pass  through  each  characteristic  point  with  slopes 


/ du\  _  -b  ±  Vb2  -  4ac 

\dr/c  =  2a 

where 


(69) 


(70) 
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0r2  2r2  +  2K2' 

c  c 


'lu2  p[Er3  +Kw  r  "ju  - 

J  c  r-1  PL  e  c  cj  o 


20r2w2  -  20E2r4 
c  c  c 


(71) 


(r-lMrf+K2) 


[2(rP-l)r^  +  (2PY  +  P2  -  2)rcK2Ju^ 


(r-i)< 


= - s-  [(5r-2)0Er4  +  (5Y-2+p)0EK2r2  +  (Y+20Y-30)pKr 

r  +K2)  L  c  c 


2 

w 

c  c 


Y0K?w  lu2  +  S-g  j~3E2r3  +  3E2K2r3  +  (r3+r  K2-2pr3)w2 
cj  c  r2+K2  L  c  c  \  c  c  c/  c 


+  2pKEr 


5w  lu 
C  CJ  c 


(72) 


If,  for  example,  r  and  u  are  specified,  the  other  coordin- 
c  c 

ates  T  and  w  of  the  characteristic  point  are  found  from  Eqs.  (65) 
c  c 

and  (66),  respectively*  The  slopes  of  the  two  solutions  are  then 
found  from  Eq.  (69)*  This  information  is  sufficient  to  start  inte¬ 
gration  of  Eqs.  (59),  (62)  and  (63)  at  the  characteristic  point  in  both 
the  upstream  and  downstream  directions. 

In  Reference  2,  a  study  is  made  of  an  accelerator  similar 
to  the  one  shown  in  Figure  2.  The  electrodes  are  assumed  to  be  para¬ 
llel  to  the  r-0  plane  so  that  at  these  surfaces  w  =  0.  Eqs.  (59)  and 
(63)  were  solved  for  supersonic  flows  with  the  assumption  w  =  0  every¬ 
where.  It  is  claimed  this  is  a  reasonable  approximation  as  long  as 
K  <  1.  The  real  criteria  should  be  how  large  is  the  term 

r  +  n 

in  Eq.  (62).  Even  if  K  <  1  this  term  might  be  large  and  make  the 
approximation  .a  poor  one.  In  order  to  obtain  an  exact  solution  one 
must  return  to  the  set  of  Eqs.  (5l)”(5**)  which  allow  for  a  variation 


in  the  z  direction 
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IV.  CASE  2 


The  non-dimensional  forms  of  Eqa.  (38)  and  (39)  are 


d  =  -T~3  |-wr  +K(Er  +  u)  (73) 

r  r  +VL  L  J 


Jz  =  ”T72[Er2+w+Kw]  (74) 

r  +K 


The  introduction  of  w  and  the  Hall  current  usually  results  in  a 
radial  component  of  current.  If  we  specify  that 


K(Er  -i-  u) 
r 


(75) 


so  that  =  0,  Eqs.  (59) *  (62)  and  (63)  are  satisfied  by  the  solu¬ 
tion 


u 


-  NE 
N  +  m 


r 


(76) 


w 


P 


KE 

P  +1 

-  E2m 
P+1 


(77) 

(78) 


m(P  +  1) 
2 

pEr 


The  Mach  number  is  constant 

«  -  (t4t> 


(79) 


(80) 


This  particular  solution  reduces  to  Case  C  when  K  a  0 
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V.  CASE  3 

If  K  =  0,  Eqs.  (59M63)  reduce  to 


du 

dr 


2y_  2  4y-2  0-2  2 

"  -4-  u  —  ■ 1  ~s  1  Eru  -  2E  r 
V-l  Y-l 


P u2‘T) 


F 

G 


(81) 


w  = 


:lr 


-P 


(82) 


4^  a  2pE[Er  +u]  -  2u  — ■  -  2w  ~  (83) 
or  dr  dr 

where  ^  is  a  constant  of  integration.  With  the  Hall  current  absent, 
solutions  can  still  be  obtained  with  /  0.  The  numerical  integra¬ 
tion  of  these  equations  is  treated  in  the  same  way  as  in  Case  1. 
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VI.  CASE  4 


If  =  0,  so  that  w  =  0,  Eqs.  (8l)  and  (83)  become 


du  _  21 

[  2r  2 

— T-  U  + 

-Y-l 

4y-2  ^  „„2  2 

— j-  Eru  +  2E  r 

+  uT 

dr 

r 

r  2  2 

u  -  T 

Y-l 

(84) 


dT 

dr 


2|3E(Er  +u)  -  2u 


(85) 


In  this  case,  the  flow  velocity  is  in  the  radial  direction  and  the 
electrodes  shown  in  Figure  2  are  planes  perpendicular  to  the  wire*  It 
is  convenient  to  split  this  problem  into  two  subcases,  depending  on 
the  sign  of  E. 

If  E  <  0,  the  introduction  of  the  similarity  parameter* 

C  =  -  Er  (86) 


dT 

dr 


2p  (  C  -  u)  -  2u  £ 


(88) 


Again  we  restrict  our  interest  to  solutions  that  do  not  terminate  in 


•  **u  T 

The  similarity  parameters  U  =  -=-  and  t  =  — w  could  also  be  used.  These 

parameters  are  used  in  the  treatment  of  the  constant  area  channel  in 
Appendix  A. 


-  16  - 


limit  lines.  Since  w  =  0,  the  flow  is  sonic  at  the  limit  circle,  so 
the  solutions  we  seek  must  pass  through  the  sonic  velocity. 

The  singularity  condition  G  =  0  is  the  same  as  Eq.  (65) » 
while  the  regularity  condition  F  =  0  leads  to 


1  -2y 


(89) 


The  two  roots  ^  and  are  plotted  as  a  function  of  p  in  Figure  6  for 

a  monatomic  gas.  Only  positive  values  of  t  are  plotted  since  u  must 

c 

be  positive.  For  this  reason  is  not  shown  for  p  <0.6*  The  graph 
shows  that  0  s  i|f^  ^  0.4  and  1.0  <  . 

Eq.  (89)  shows  the  locus  of  the  projection  of  the  character¬ 
istic  points  (C  ,  u  ,.T  )  on  the  u  -C  plane  is  a  straight  line  for 
c  c  c  c  c 

each  value  of  p.  Some  of  these  lines  are  shown  in  Figure  7-  No  charac 
teristic  points  lie  in  the  region  separated  by  the  lines  p  =  ». 

Following  the  method  used  in  Case  1,  the  slopes  of  the  two 
solutions  passing  through  a  characteristic  point  are 


-  b  ± 


2a 


where 

a  =  8  If 

b  =  (I5p -lH2 -12pt|i 
c  =  (2-5t)pt 

For  each  value  of  two  values  of  (du/dC) 

c 


(90) 


(91) 


are  obtained.  The  slopes 
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are  a  function  of  (3  only,  and  hence  are  constant  along  each  of  the 
straight  lines  in  Figure  7.  The  variation  of  slope  with  (3  is  shown 
in  Figure  8*  Only  one  of  the  four  curves  is  negative,  so  that  most 
of  the  solutions  are  accelerating  flows*  The  curves  corresponding  to 
+2  terminate  at  f3  =0.6  where  is  infinite,  and  the  curves  corres¬ 
ponding  to  terminate  at  p  =0.3  where  the  slopes  become  complex. 
These  results  are  shown  in  Figure  7* 

Numerical  solutions  were  obtained  for  C  =  1.0  and  (3  =  2.0. 

c 

The  results  are  compared  with  the  ordinary  gasdynamic  flow  E  =  B  =  0 
in  Table  III. 


TABLE  III 


The  solid  curve  in  Figures  9-15  corresponds  to  an  accel¬ 


erating  flow  which  becomes  asymptotic  to  the  solution  given  by 
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Eqs.  (76)  -  (80)  in  Case  2  with  K  =  0.  The  solution  begins  at  a  finite 
radius  with  infinite  density,  finite  pressure,  and  zero  velocity,  tem¬ 
perature  and  Mach  number.  Energy  is  continuously  fed  into  the  gas  as 
it  accelerates  through  the  sonic  velocity. 

The  solution  shown  by  the  dashed  curves  is  an  accelerating 
flow  whose  temperature  and  pressure  increases  as  the  Mach  number  and 
density  decreases.  The  solution  begins  with  zero  pressure  and  density,  . 
and  finite  velocity,  temperature  and  Mach  number-  As  the  velocity 
increases,  the  rapid  rise  in  temperature  causes  a  drop  in  the  Mach 
number  until  the  initially  supersonic  flow  becomes  subsonic.  Again 
energy  is  continuously  added  to  the  gas,  but  at  any  point  in  the  flow 
this  solution  has  a  higher  stagnation  enthalpy  than  the  previous  one. 

The  dot-dashed  curves  begin  at  C  =  0  with  infinite  pressure 
and  density,  and  finite  velocity,  temperature  and  Mach  number.  The 
flow  accelerates  through  the  sonic  velocity  and  becomes  asymptotic  to 
the  solution  for  Case  2.  In  order  to  accomplish  this,  a  minimum  occurs 
in  the  temperature  curve  and  a  maximum  in  the  Mach  number  curve.  Energy 
is  extracted  from  the  flow  up  to  £  ='  2.2  and  added  to  the  flow  beyond 
that  point.  Thus,  the  flow  acts  first  as  a  generator  and  then  later 
as  a  motor. 

The  solution  plotted  as  heavy  dashed  curves  is  the  only  de¬ 
celerating  flow  at  the  sonic  point.  The  flow  begins  at  a  finite  radius 
with  zero  temperature  and  pressure,  and  finite  velocity,  density  and 
Mach  number.  The  supersonic  flow  decelerates  through  the  sonic  velo¬ 
city,  reaches  a  minimum  velocity,  and  begins  to  accelerate.  The  density 
reaches  a  maximum  while  the  flow  is  still  decelerating.  The  flow 
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changes  from  a  generator  to  a  motor  after  passing  the  sonic  point* 

The  second  subcase  occurs  when  E  >  0.  Introducing  the  simil¬ 
arity  parameter 


£  =  Er 


Eqs*  (84)  and  (85)  become 


(92) 


§  -  2e«  +  u)  -  2u  If 

Following  the  procedure  used  in  the  first  subcase, 


where 

a  =  89 

b  =  (15P  -  1)<P2  +  12ptp 

c  =  (2+5?)?? 


(93) 


(94) 


(95) 


(96) 


The  variation  of  the  quantity  9  appearing  in  Eq.  (94)  is  shown  in 
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Figure  6«  This  quantity  is  positive  only  if  (3  <0.6,  and  takes  on  all 
positive  values.  Hence,  all  points  in  the  u-£  plane  can  be  charac¬ 
teristic  points.  The  slopes  of  the  solutions  at  the  characteristic 
points  are  plotted  in  Figure  14  as  a  function  of  p.  The  slopes  are 
both  negative  and  are  real  for  0«32  <  (3  <  0«60.  The  locus  of  possible 
characteristic  points  which  yield  real  solutions  are  shown  in  Figure  15* 

Numerical  solutions  were  obtained  for  £  -  1.0  and  |3  a  0.5* 

c 

These  results  appear  in  Figures  16-20.  Only  curve  I  is  shown  because 
an  infinite  number  of  curves  pass  through  the  characteristic  point  with 
the  same  slope  as  curve  II,  which  is  the  envelope  of  the  curves.  The 
numerical  integration  in  a  direction  away  from  the  characteristic  point 
was  found  to  be  unstable  while  integration  towards  the  characteristic 
point  was  stable.  Two  solutions  called  Ila  and  lib  are  shown  in  Fig¬ 
ures  21-25  for  a  small  region  near  the  characteristic  point.  Away  from 
the  characteristic  point  these  solutions  are  practically  the  same  as 
curve  I  and  cannot  be  distinguished  in  the  scales  used  in  Figures  16- 
19*  No  attempt  was  made  to  extend  the  envelope  II  or  to  extend  the 
curves  Ila  and  lib  into  the  region  where  £  <  1. 

Curve  I  begins  at  a  finite  radius  with  zero  temperature  and 
pressure,  and  infinite  Mach  number.  The  velocity  and  density  are 
finite.  As  the  flow  decelerates  through  the  sonic  velocity  the  tem¬ 
perature  and  pressure  increases  as  the  density  decreases.  Energy  is 
added  to  the  flow  for  all  solutions  with  E  >  0. 

The  reason  why  an  infinite  number  of  solutions  must  pass 

through  the  characteristic  point  is  that  both  curves  I  and  II  are 

decelerating  and  supersonic  when  £  <  £  ,  and  both  decelerating  and 

c 


2i 


subsonic  when  £  >  5  «  This  means  they  pass  from  region  1  to  region  3 
in  Figure  26*  The  curve  AC  is  the  locus  of  characteristic  points,  or 
the  intersection  of  the  surface  F  =  0  with  the  sonic  surface  G  =  0. 

In  regions  1  and  2  above  the  sonic  surface  the  flow  is  supersonic,  and 

in  regions  3  and  below  the  sonic  surface  the  flow  is  subsonic.  The 

curves  I  and  II  intersect  at  the  characteristic  point  B.  The  surface 
E  is  generated  by  curve  I  as  B  moves  from  A  to  C  and  the  surface  D  is 
a  similar  surface  generated  by  curve  II#  Solutions  which  pass  through 
points  in  region  1  lying  between  surfaces  D  and  E  must  pass  through 
the  curve  AG  with  the  two  permissible  slopes  of  curves  I  and  II#  The 
solutions  cannot  penetrate  the  surfaces  D  or  E  except  along  AC,  for 
this  would  produce  a  characteristic  point  which  can  occur  only  along 
AC.  Neither  can  they  be  parallel  to  AC  because  the  projection  A'B' 
on  the  u-£  plane  has  a  positive  slope  while  ^  <  0  in  region  1.  Thus, 

as  numerical  integration  shows,  curve  II  is  an  envelope  of  solutions 

such  as  Ila  and  lib  which  pass  through  the  characteristic  point  with 
the  same  slope. 

In  all  the  cases  for  E  <  9,  which  are  shown  in  Figure  9, 
the  two  curves  passing  through  the  characteristic  point  are  always 
separated  by  a  sonic  surface  and  the  surface  F  =  0  for  any  given  value 
of  £.  Thus,  at  the  same  £  one  solution  is  subsonic  while  the  other 
is  supersonic#  This  means  the  surfaces  corresponding  to  D  and  E  in 
Figure  26  always  lie  in  different  regions,  and  only  two  curves  can 
pass  through  the  characteristic  point# 

Some  solutions  for  accelerating  supersonic  flows  were  pre¬ 
sented  in  Reference  3*  One  type  of  solution  is  asymptotic  to  Case  2 
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with  K  =  0.  The  second  type  of  solution  accelerates  to  a  maximum 
velocity  and  then  decelerates  and  terminates  at  a  limit  line*  No 
solutions  for  acceleration  through  the  sonic  velocity  were  presented. 
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VII.  CASE  5 


This  case  deals  with  solutions  for  finite  values  of  R  . 

m 

First,  the  method  of  small  perturbations  is  used  to  find  solutions 

which  hold  for  small  values  of  B  .  The  one  solution  of  Case  A  which 

m 

passes  through  the  sonic  velocity  is  chosen  as  the  zero  order  solution 
Let 


u 

P 

P 

B0 


R  u. 

2 

+  R  u_  +  . 

m  1 

m  2 

_2 

R  P-. 
nr  1 

+  R  P-  +  • 
nr  2 

R  P, 
m  l 

^2 

+  R  p_  +  . 
m  d. 

*  sJih\  * 

With  E  =  —  =  0,  Eqs •  (42)-(50)  reduce  to 


rpu  =  m 

(98) 

du  dp  XT 

pu  d7  +  a  d7  =  "  NuBe 

(99) 

2 

p  +  pu  .  =  p 

(100) 

dBfi  Bfi 

— —  +  —  =  R  uB„ 

dr  r  mo 

(101) 

Substitution  of  Eq.  (97)  into  Eqs.  (98)-(101)  and  collection  of  like 
powers  of  Rm  yields: 


rp  u  =  m 
o  o 


(102) 


p  u 


(103) 
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p  +  p  u  =  p 

*0  o  o  o 

<Be)  (B«) 


dr 


o  o 

+  — - 


a  0 


(£> 


pluo  +  poul  =  0 


dui  r  ^du0  dpi 

poUo  "dT  +  Wl  +  PlUoJ  “dr“  +  a  "dr"  = 


-N 


2 

<0  *2“o(B8)(Be). 


O  '1J 


pl  +  2p0U0Ul  +  P1U0  =  P1 


<0,  CO,  ,  ^ 

sr-'-r-  -  uoCv. 


(104) 


(105) 


(106) 


(107) 

(108) 


(109) 


In  this  analysis  the  higher  order  terms  will  be  neglected. 

The  zero  order  equations  (102)-(106)  are  satisfied  by  the 
solution  of  Case  A  with  p  =  y  •  . 


(Be)„  ■  SS 

P°  ■  (}-$ 


(no) 


/ 
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Combining  Eqs.  (106)  and  (108), 

pi  ■  -  (b  *  £>  ui 


(112) 


and  integration  of  Eq.  (109)  leads  to 

CO 


in  C1R 


(113) 


where  is  the  integration  constant.  Elimination  of  and  (j 3^) 
from  Eq.  (107)  yields 


dOhij^) 

dR 


-  4in  C  R 


_  F 


r[(y+1)-(y-1)E2]  g 


(114) 


The  singularity  condition  G  =  0  is 

-  m 

and  the  Regularity  condition  F  -  0  is 

C1  =  45 

Integration  of  Eq.  (114)  results  in 

.45  R 


(115) 


(116) 


-4 

(Y+1)h 


/ 

45 


in  x 


x(l-x2) 


dx 


(117) 


The 


variation  of  CO’  v  "i  and  p1  with  R  when  y  =  5/3  is 
shown  in  Figure  27.  (j3g^  is  negative  upstream  of  the  sonic  point,  and 
zero  at  the  sonic  point.  Downstream  of  the  sonic  point  it  is  positive, 
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and  reaches  a  maximum  at  R  =  2e.  The  curve  for  u^  increases  rapidly 
from  zero  to  a  maximum  at  the  sonic  point.  The  perturbation  in  pressure 
and  the  perturbation  in  density  decreases  from  zero  to  -  2»5» 

The  previous  method  is  valid  only  for  small  values  of  the 
magnetic  Reynolds  number.  The  application  of  the  method  to  some  other 
solutions  of  Case  A  other  than  the  one  with  (3  =  y  ^  is  difficult  because 
of  the  complicated  form  of  the  zero  order  solution.  For  these  reasons 
an  attempt  is  made  to  solve  Eqs.  (98)-(lOl)  directly  for  finite  values 
of  R^.  The  numerical  approach  of  Case  4  is  used. 

Integration  of  Eq.  (101)  yields 


(118) 


As  will  be  shown  later,  all  solutions  passing  through  the  sonic  velo¬ 
city  are  decelerating  so  that  r^  is  either  zero  or  the  place  where 
u  =  1.  Elimination  of  the  density  and  pressure  from  the  momentum  equa¬ 
tion  results  in  the  expression 


du 

dr 


-PrVB2  ,g)  u(l-u2) 
r  7^2^1777 


(119) 


Restricting  our  attention  to  y  =  5/3,  and  substituting  Eq.  (118)  for 

r 

V 


-5u^(3  exp  2R  f  u  dr  +  u(l-u2) 

du  _ L  m  ri  J  =  F 

dr  =  r[4u2-l]  '  ® 


(120) 


Since  we  wish  to  start  integration  of  Eq.  (120)  at  the  sonic  point,  it 
might  appear  that  the  unknown  integral  in  the  numerator  might  be 
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troublesome*  Fortunately,  this  is  not  the  case  because  the  regular¬ 
ity  condition  requires  the  product  of  p  and  the  exponential  factor  to 
be  a  known  constant  at  the  sonic  point.  Thus,  it  is  possible  to  inte¬ 
grate  without  knowing  the  value  of  Bq  at  the  sonic  point.  Bq  is  cal¬ 
culated  from  Eq.  (118)  after  u  is  found.  The  singularity  condition  is 

u  s  £  and  the  regularity  condition  leads  to 
c 


0 


=3 


1 

5 


exp 


u  dr 


(121) 


so  that  Eq.  (120)  becomes 


-|l  +  3  exp 

rT 

2R  /  u  dr 

m  1 

y + u 

du  V 

Ui  f 

Jr 

-  c  -J 

/ 

dr 

r 

(4u2  - 1) 

(122) 


The  sonic  point  is  now  the  logical  point  to  begin  numerical  integration. 
The  slope  at  the  sonic  point  is  given  by 

/,  \  .  -1  ±  Jl  -  R  r  ' 

fe)  -  — V  mc  (123) 

\  'c  c 

The  slope  is  always  negative  and  real  if  Rmrc  s  2/3 *  This  equation  is 

plotted  in  Figure  28  for  several  values  of  R  .  Through  each  character- 

m 

istic  point  (0.5,  r  )  there  can  be  at  most  two  possible  slopes,  both  of 

c 

which  are  negative.  The  larger  of  the  two  slopes  (in  absolute  value) 
lies  in  Region  I.  Each  slope  in  Region  II  corresponds  to  an  infinite 
number  of  solutions  as  in  Case  5  when  E  >  0.  This  was  expected  since 
both  slopes  are  negative  and  are  also  subsonic  and  supersonic  in  the 
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same  regions* 

If  we  limit  our  calculations  to  slopes  in  Region  I,  the 

solutions  which  pass  through  rc  =  1  appear  like  those  in  Figure  29* 

As  R  is  increased  from  0  to  the  maximum  value  of  2/3,  the  value  of  3 
m 

drops  from  0.6  to  0.31*  All  the  curves  choke  at  a  limit  circle  whose 

radius  decreases  rapidly,  as  R^  -*  2/3« 

By  fixing  R^  and  integrating  backwards  from  each  value  of 

rc ,  r^  and  p  are  obtained.  These  results  are  plotted  in  Figure  30  for 

slopes  in  Region  I.  For  R^  =  0,  p  =  0.6*  As  R^  increases,  p  can  take 

on  values  less  than  0.6.  The  maximum  value  of  r^  decreases  with  an 

increase  in  R  . 

m 

The  line  r.  =  3*5  intersects  the  curve  R  =  0.5  in  points  A 
i  m 

and  B.  These  points  also  appear  in'  Figure  28  in  Region  I.  Note  that 
point  A  does  not  lie  on  the  dashed  line,  which  divides  the  two  regions, 
so  that  a  portion  of  the  curve  extends  below  point  A  in  Figure  30. 

Points  between  A  and  B  in  Figure  30  have  values  of  p  which  correspond  to 

solutions  which  pass  through  the  sonic  velocity  with  slopes  in  Region  II. 

Therefore,  as  we  increase  p,  there  is  first  a  jump  in  slope  from  A  to  A', 
in  Figure  28,  followed  by  a  path  to  B'  along  the  Rm  =  0.5  curve,  which 
is  terminated  by  a  second  jump  to  B.  If  r^  is  reduced  to  2.5,  point  D 
lies  in  Region  I  and  all  points  from  D  to  C,  including  C,  lie  in 
Region  II.  Points  which  do  not  lie  under  the  curve  Rm  =  0.5  choke. 

In  Figure  31,  the  range  of  values  of  Rm  and  p  which  permit 

passage  through  the  sonic  velocity  is  shown  for  r^  =  0.35*  The  solid 

curve  CED  contains  points  in  Region  I,  and  the  shaded  area,  including 
the  dashed  curve  CD,  contains  points  in  Region  II.  Line  segment  AB  is 


-  29  - 


the  same  as  in  Figure  30. 

The  solution  for  r^  =  0.35  is  shown  in  Figure  32.  As  0  is 
reduced  to  the  value  0.41  all  the  solutions  choke.  For  p  =  0.41,  we 
have  point  B  in  Figures  30  and  31  and  is  in  Region  I  of  Figure 

28.  As  p  is  decreased  the  slope  jumps  to  Region  II  from  B'  to  A'  , 
and  then  back  to  Region  I  for  point  A.  For  values  of  (3  <  0.31  which 
corresponds  to  point  A  the  solutions  again  choke.  An  infinite  number 
of  solutions  pass  through  the  sonic  velocity  in  the  range  0.31  s  P  s  0.41. 

If  Rm  =  0,  this  region  shrinks  into  point  D  in  Figure  31 
and  p  =  0-6*  The  solution  is  shown  in  Figure  33*  A  comparison  can  be 
made  with  Figure  3  where  the  stagnation  boundary  condition  was  used 
instead  of  r^  =  const. 


VIII.  CASE  6 


Case  A  treats  the  limiting  case  =  0,  N  jl  0,  in  which  the 
magnetic  field  interacts  with  the  flow,  but  is  not  itself  changed  by 
the  flow.  For  non-zero  values  of  Rm  and  N  both  the  flow  and  magnetic 
fields  mutually  interact.  This  case  is  concerned  with  the  other  limit¬ 
ing  case  when  N  =  0,  R^  /  0.  In  this  case,  the  flow  distorts  the 
magnetic  field  but  is  not  itself  changed  by  the  interaction. 

Assuming  r~  =  0,  Eqs.  (42)  -(48)  reduce  to 


rpu  =  m 


(124) 


du  _  dp 

pu  d7  +  a  d5  = 


(125) 


H  =  1 


(126) 


d_B  B 
dr  r 


R  (E  +  uB) 
m . 


(127) 


The  first  three  equations  are  the  ordinary  gasdynamic  equations 
which  integrate  to  yield 


r 


m 


u(l-u2)5/2 


(128) 


This  is  the  source  flow  which  corresponds  to  the  ^  =  0  curve  in  Fig¬ 
ure  3* 


Equation  (127)  can  be  placed  in  the  form 
ER 


u  .  T  r 

R  f  r  -r—  exp  -R  f  udr 
mJ  du  mj 

u  r 


du  +  B  r  ■ 
c  c 


B 


(129) 
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where  is  given  by  Eq.  (128)  with  =  0«5*  Substitution  of  Eq.  (128) 
into  Eq.  (129)  yields 


B 


ER  J  f(u)g(u,mR  )du  + 1 
m*l/2 _ m 

~  g(u,mR  ) 
m  m 


(130) 


In  Eq.  (130)  the  boundary  condition  B  r  =  m  has  been  applied. 

c  c 

The  results  of  numerical  integration  of  Eq.  (130)  are  plotted 

in  Figures  (34)  -  (37).  For  a  subsonic  source  with  E  =  0,  as  r/m  -*  00 , 

B  -*  0  when  mR  <  1,  B  -*  O.366  when  mR  =1,  and  B  -•  »  when  mR  >1. 

These  results  are  shown  in  Figure  (34).  If  E  is  increased  to  one, 

B  -*  0  when  mR  =  0,  and  B  const,  when  0  <  mR  s  1.  When  mR  >  1, 

.  m  m  m 

B  -»  «.  These  results  appear  in  Figure  (35)* 

For  a  supersonic  source  the  magnetic  field  B  *•  »  with  r/m 

when  mR  >0  and  B  -*  0  when  mR  =  0.  The  distribution  of  B  is  shown  in 
m  m 

Figure  36  for  E  =  0  and  in  Figure  37  for  E  =  1.0. 

The  limiting  values  of  B  for  large  r/m  are  summarized  in 


Table  IV. 


33 


IX.  CONCLUSIONS 

Solutions  have  been  obtained  for  magnetogasdynamic  source 
flow  with  crossed  electric  and  magnetic  fields.  The  solutions  have 
been  restricted  to  those  which  pass  through  the  sonic  velocity. 

For  the  case  where  Rffl  =  K  =  0,  the  sonic  circle  can  be  placed 
at  any  radius.  If  £  >  -0,  an  infinite  number  of  solutions  decelerate 
through  M  =  1  with  the  same  velocity  gradient.  There  is  one  other  solu¬ 
tion  which  also  decelerates  through  M  =  1  with  a  different  gradient  in 
velocity.  If  E  <  0,  four  different  solutions  can  pass  through  sonic 
velocity  at  any  given  location  of  the  sonic  circle.  Three  of  these  are 
accelerating  flows. 

The  effect  of  non-zero  magnetic  Reynolds  number  was  found  for 

source  flow  with  E  =  0.  When  R^  =  0  there  is  one  solution  which  passes 

through  M  =  1  when  p  =  0.6*  When  R^  >  0  there  is  a  range  of  p  in  which 

solutions  are  possible.  If  R  is  increased  beyond  a  certain  limit,  no 

m 

solutions  are  possible. 

The  displacement  of  the  magnetic  field  by  ordinary  gasdynamic 

source  flow  was  obtained  for  the  limiting  case  when  p  =  0,  R  ^  0.  For 

m 

a  supersonic  source  the  magnetic  field  always  increases  to  infinity 

with  the  radius  for  finite  R  .  For  a  subsonic  source  this  occurs  only 

in 

if  mR  >  1.  For  other  values  of  mR  the  field  approaches  a  constant 
m  m  ** 

value  or  zero. 

These  results  indicate  many  possibilities  exist  for  the 
passage  through  the  sonic  velocity  without  the  need  of  a  throat  in 
magnetogas dynamics. 
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APPENDIX  A 

ACCELERATION  THROUGH  A  CONSTANT  AREA  CHANNEL 

This  section  deals  with  the  flow  through  a  constant  area 
channel  with  crossed  electric  and  magnetic  fields*  The  configuration 
is  shown  in  Figure  A-l* 

The  governing  equations  in  dimensionless  form  are 

pu  =  m  (A-l) • 

-N(E+u)  (A-2) 

»  2NE(E  +  u)  (A-3) 


du  dp 

pu  rr  +  a  *“• 
dx  dx 


fdT  .  du  1 

‘’''te  * 2u  sj 


pT 


Eliminating  the  pressure  and  density  yields 


where 


dU 

dx 


-BP(1  -P)(2  -5U) 


3U2  -  t 


£  =  2P(1-U)  -  2yg 


F 

G 


The  singularity  condition  G  »  0  leads  to 


and  the  regularity  condition  F  =  0  leads  to 


U 

c 


=  0, 


D 

c 


=  1, 


U  » 

c 


(A-5) 

(A-6) 


(A-7) 


(A-8) 


0.4 


( A*9) 


The  two  slopes  at  a  characteristic  point  are 


(A. 10) 


The  properties  at  the  three  characteristic  points  are  summarized  in 
Table  A-I. 


TABLE  A- I 


u 

c 

t 

c 

( M ) 

C 

0 

0 

0 

1.5P 

20 

20 

0*4 

0.48 

0 

0.75P 

1.20 

0.60 

1.0 

3-0 

0 

-0.375P 

0 

0.750 

Since  Eqs.  (A-5)  and  (A-6)  do  not  contain  x,  we  can  divide 
the  two  equations  to  eliminate  dx  to  obtain 


dU 

dt 


P(5P-2) 

-  10U3  +  10U2  -  2t 


(A. 11) 


This  equation  in  the  phase  plane  is  independent  of  0.  Two  character¬ 
istic  points  exist  and  the  properties  of  the  solutions  through  these 
points  are  summarized  in  Table  A-II. 


TABLE  A-II 


u 

c 

*c 

0 

0.4 

0 

0.48 

0 

0 

0.75 

1.25 

-  36 


The  solutions  in  the  phase  plane  are  shown  in  Figure  A- 2. 

The  projection  of  the  solutions  on  the  U-x  and  t-x  planes  are  shown 
in  Figures  A-3  and  A-4,  respectively. 

Solution  la  accelerates  from  zero  temperature  and  velocity 
through  the  sonic  velocity  to  M  *  at  x  =  "•  The  temperature  reaches 
a  maximum  at  x  *  0«9  where  U  =  0.9*  Energy  is  always  added  to  the 
flow,  but  the  rate  of  addition  diminishes  with  x.  In  solution  lb,  the 
velocity  is  constant  and  the  temperature  increases  linearly  with  x, 
so  that  the  initially  supersonic  flow  becomes  subsonic.  Energy  is 
added  to  the  flow  at  a  constant  rate. 

Energy  is  first  extracted  from  the  flow  when  it  is  super¬ 
sonic  and  then  added  when  it  becomes  subsonic  in  solution  Ha.  The 
flow  begins  at  U  =  2.34  with  zero  temperature  and  decelerates  to 
U  »  0.4  and  infinite  temperature.  Solution  lib  is  a  constant  Mach 
number  flow  at  the  sonic  velocity. 

A  closed  form  solution  corresponding  to  la  and  lb  appears 
in  Reference  4.  Solutions  la  and  Ila  were  discussed  in  Reference  5, 
but  calculations  were  not  made* 
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Fig.  1.  Coordinate  system  and  geometry  of  MGD  source. 


Fig.  2 


Coordinate  system  and  geometry  of  diverging  channel 


Case  A.  Velocity  distribution  for  various  values  of 
the  interaction  parameter  (3.  Isoenergetic  flow  E  =  0 
P  >  0,  source;  f3  <  0,  sink. 


Fig.  5a.  Location  of  sonic  circle  and  limit  circle  when  solution 
is  inside  the  limit  circle. 


Fig.  5b.  Velocity  distribution  for  solution  inside  limit  circle. 


Fig.  5c.  Location  of  sonic  circle  and  limit  circle  when  solution 
is  outside  the  limit  circle. 


Fig.  5d.  Velocity  distribution  for  solution  outside  limit  circle 


Fig.  6 


Variation  of  1^,  t2  and  0  with  (3 


Case  4,  E  <  0 
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Fig.  7*  Locus  of 
Case  4, 


<  0 
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Fig.  10.  Variation  of  T  with  £  for  £  =1  and  j3  =  2.  Case  4 
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Fig.  13.  Variation  of  ^  with  C  for  C  =  1  and  B  =  2.  Case 

me 
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Fig.  16*  Variation  of  u  with  £  for  £  =1  and  0  =  0.5»  Case  4, 
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E  >  0. 


Fig.  21. 
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E  >  0 
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Fig.  23.  Variation  of  M  with  £  for  £  =  1  .nd  P  =  0.5*  Case  4, 
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E  >  0 


3hs 
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Fig.  24.  Variation  of  £  with  r  for  r  =  1  and  (3  =  0.5*  Case  4, 


E  >  0 


Q.IE 


solution  in  the  u  -  T  -  £  coordinate  system 


Fig.  27.  Variation  of  the  perturbation  variables  with  R. 


Case  5 


29*  Velocity  distribution 
of  Figure  28*  Case  5 


from  region 


Variation  of  r 


Fig.  Jl.  R  vs.  0  when  r.  =  0.35* 
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Fig.  36*  Magnetic  field  distribute 
=  0. 


E 


for  supersonic  source  with 


Fig.  37.  Magnetic  field  distribution 


Fig.  A-2.  Solution  of  constant 
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Fig.  A-3*  Distribution  of  the  temperature  function  t 


